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ABSTRACT 

The cores of Arp 220, the closest ultraluminous infrared starburst galaxy, provide 
an opportunity to study interactions of cosmic rays under extreme conditions. In this 
paper, we model the populations of cosmic rays produced by supernovae in the central 
molecular zones of both starburst nuclei. We find that ~ 65 — 100% of cosmic rays 
are absorbed in these regions due to their huge molecular gas contents, and thus, the 
nuclei of Arp 220 nearly complete proton calorimeters. As the cosmic ray protons 
collide with the interstellar medium, they produce secondary electrons that are also 
contained within the system and radiate synchrotron emission. Using results from 
X 2 tests between the model and the observed radio spectral energy distribution, we 
predict the emergent 7-ray and high-energy neutrino spectra and find the magnetic 
field to be at milligauss levels. Because of the extremely intense far-infrared radiation 
fields, the 7-ray spectrum steepens significantly at TeV energies due to 7-7 absorption. 

Key words: neutrinos - cosmic rays - galaxies: individual: Arp 220 - galaxies: 
starburst - gamma rays: galaxies - radio continuum: galaxies 


1 INTRODUCTION 


Arp 220 is the nearest (d = 77 Mpc) example of an ultralu¬ 
minous infrared galaxy (ULIRG) that supports star forma¬ 
tion at extreme levels. It contains two nuclei separated by 
350 pc, both surrounded by mas sive discs of dense molecu- 

lUlIiijL 


lar gas (e.g., 


Mundell et al 


Smith et al. 


Downes fe Solomon! 19981: 


200 ll : Sakamo to et al.l 20081 : Scoville et al.l 


2015|). Radio detections of supernovae at a rate of 1-3 yr 


(jSmi^e^alJ I 19981 : iRovilos et al.ll2005l : iLonsdale et al.ll2006l : 
iParra et alj 20071 1 confirm that huge populations of massive 
stars are present with an implied star formation rate (SFR) 
of > 100 Mq yr _1 . Although Arp 220 could contain active 
galactic nuclei (AGNs), particularly in the western nucleus, 
the observed supernova rates indicate that star formation 
provides a substantial fraction of the power radiated by the 
nuclei. 

The nuclei of Arp 220 provide access to the high- 
intensity mode of star formation in dense molecular me¬ 
dia that appears to have been more common in young 
galaxies. These types of environments are of special inter¬ 
est from a range of perspectives, including the information 
they can provide regarding the role of galactic winds, cos¬ 
mic rays, and magnetic fields in feedback processes that in- 
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fluence galaxy evolution. Previous investigations show that 
Arp 220 is likely to be a hadronic cosmic ray calorimeter 
where all of the power in cosmic rays is abso rbed within 
the nuclear starburst zones (e.g., iTorresl [20041 1. Both nu¬ 
clei also contai n extremely intense far-infrared (FIR) radia¬ 
tion fields (e.g ., Soifer et al . 200d : [Barc.os-Munoz et al.ll2015l : 
IScoville et al.ll2015l ). and the west nucleus i s optically thick 
in th e FIR to wavelengths of A « 800 fi dSakamoto et al .1 
120081 : IPapadopouIos et ahl 1201(f ). The production of the 
observed radio synchrotron emission then requires mag¬ 
netic fields of mil ligauss strength (e.g.. iLisenfeld et al .11 1 99(il : 
lLacki et al.ll2Qloh . 


In this paper, we study cosmic ray interactions in the 
Arp 2 20 starburst nuclear reg ions using an updated version 
of the I Yoast-Hull et al.l (12013 ) models, hereafter YEGZ. We 
develop a model with two spatial zones to accurately repre¬ 
sent the inner and outer regions of the western nucleus a s 
defined by its molecular gas properties dAalto et al.|[2009l ). 
We incorporate photopion energy losses and photon-photon 
interactions to account for the extreme FIR radiation field. 
We calculate the hadronic calorimetry fraction for each nu¬ 
cleus for the best-fitting radio models, and we predict the 
total 7 -ray and neutrino fluxes. 

In Section 2, we review the physical parameters which 
we selected for the models. Section 3 details the basic as¬ 
sumptions of the models and our findings for the Arp 220 
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2 Yoast-Hull, Gallagher, & Zweibel 


Table 1. Input Model Parameters 


Physical Parameters 

East Nucleus 

West ST 

West CND 

References 

Distance 

77.0 Mpc 

77.0 Mpc 

77.0 Mpc 


CMZ Radius 

70 pc 

90 pc 

30 pc 

1,2,3 

CMZ Disc Scale Height“ 

40 pc 

40 pc 

40 pc 

4 

Molecular Gas Mass 

6 x 10 8 M© 

4 x 10 8 M© 

6 x 10 s M© 

2,5 

Ionized Gas Mass b 

3 x 10 6 M© 

2 x 10® M© 

3 x 10® M© 


Average ISM Density“ 

~7700 cm -3 

~3500 cm -3 

~42 000 cm -3 


FIR Luminosity 

3 x 10 11 L© 

3 x 10 11 L© 

6 x 10 11 Lq 

2 

FIR Radiation Field Energy Density 1 * 

40 000 eV cm -3 

27 000 eV cm _; 

1 440 000 eV cm" 3 


Dust Temperature 

90 K 

50 K 

170 K 

2,6 

SN Explosion Rate (z'sn) 

0.7 yr- 1 

0.7 yr" 1 

1.3 yr' 1 

7 

Star Formation Rate (SFR) d 

65 M© yr' 1 

65 M© yr" 1 

120 M© yr” 1 


SN Explosion Energy 6 

10 51 erg 

10 51 erg 

10 51 erg 


SN Energy in Cosmic Ray Protons 6 

5 - 20% 

5 - 20% 

5 - 20% 


Ratio of Primary Protons to Electrons (N p /N e ) 

50 

50 

50 


Slope of Primary Cosmic Ray Source Function 

2.1 - 2.3 

2.1 - 2.3 

2.1 - 2.3 



References — (1) ^ownes_& i ^ckart (2007); (2) ^jakamoto_e^d ii (2008); (3) ^alto^etjA (2009); (4) - Scoville ii et_al ii (2015); 
(5) _Downes_^ i So]nmon (1998); (6) i Wilsoim3t_aL (2014); (7) ^onsdale^et^aL (2006). 

“Conservative estimate; ^icovine^jtjiL (2015) find thinner, denser discs. 

8 Assumes the ionized gas mass is roughly 0.5% of the molecular gas mass to keep the volume filling fraction reasonably 
less than 100%. 

“The average ISM number density is the sum of the nuclei number densities of molecular and ionized gas, such that 
?2 ism = (2 x rn p x // IIIO ] X V) + M\ on / (m p X //©,, x V). where /r mo i and are the mean molecular weights. 

^Derived from above parameters. 

“Excludes neutrino energy. 


starburst nuclei. We present concluding remarks in Section 

4. 


2 ARP 220 PHYSICAL PROPERTIES 

Due to its extreme properties, Arp 220 has been exten¬ 
sively studied across the electromagnetic spectrum. The nu¬ 
clei of Arp 220 are of particular interest as they contain 
more than half of the total bolometric infrared luminos¬ 
ity o f the galaxy (1 — 2 x 10 12 L©; e.g., iDownes fc Eckartl 
2007; Sak amoto et al.ll2008l . and references therein). As the 
nuclei are less than 100 pc in radius, the presence of an 
AGN or a ‘hot’ starburst is required to explain the ex- 
traordinarily large surface brightness in the western nucleus 

( [Downes fc Eckartil2007i : ISakamoto et al.ll2008l : IWilson et al.1 

l2014h : however, the existen ce of an AGN has yet to be 
definitively established ('e.g.. iTunnard et al.ll2015l ). Further, 
the submillimetre observations suggest that whether or not 
AGNs are pre sent, they are not the main heating source of 
the dust (e.g.. ISakamoto et al.ll200j) . 

Estimates of the FIR luminosities of the eastern and 
western nuclei range from ~ 10 10 to ~ 2 — 3 x 10 11 and 
from ~ 2 — 3 x 10 11 L© to ~ 10 12 L©, r espectively (e.g., 
ISakamoto et al.l 120081 : IWilson et al.l 120141 ). The range on 
these luminosities is quite large due to uncertainty in the 
true sizes, inclinations, and opacities of the nuclei and their 
associated molecular disc. To keep our adopted FIR lumi¬ 
nosity in rough agreement with the observed supernova rate, 
we assume values of 3 x 10 11 L© and 9 x 10 11 L© for the 
eastern and western nuclei (see Table 1). Assuming sim¬ 
ilar ratios between the nuclei for the supernova rate and 
molecular gas content, we adopt values of ;/sn = 0.7 yr -1 , 
A/md = 6 x 10 8 M© for the eastern nucleus and !zsn = 2.0 


sumed molecular gas masses favour conservative estimates, 
other estimates of the gas c ontent suggest the m asses are as 
high as ~ 2 — 4 x 10 9 M© (IScoville et al.l [20151 ). 

CO observations of the western nucleus imply a temper¬ 
ature gradient increasing towards the centre and indicate 
significant differ ences in the physi c al conditions betwee n 
the two nuclei ( Aalto et al.l 120091 : ITunnard et al.l 120151 ') . 
IDownes fc Eckart ( 20071 ) model the western nucleus as two 
distinct dust sources - a cooler (50 K) ring surrounding a 
hotter (170 K), dense dust core. We use this two-zone model 
for the western nucleus and have adjusted our single-zone 
model to account for the differences in temperature and den¬ 
sity between the two regions (see Section 3). For the east- 
ern nucleus, we assum e a single dust temperature of 90 K 
(ISakamoto et af]|2008l) . 


3 MODELS & RESULTS 
3.1 YEGZ Model 

Previously, we developed and tested a model for cosmic ray 
interactions in the central molecular zones (CMZs) of star¬ 
formin g and starburst galaxies (YEGZ; lYoast-Hull et all 
l2014al lbh. Our single-zone model accounts for a variety of 
energy losses via interactions with the interstellar medium 
(ISM), magnetic fields, and radiation fields and for energy- 
independent advective escape via a galactic wind (see Fig. 
1). The resulting cosmic ray energy spectrum depends on 
both the total cosmic ray lifetime and a power-law injec¬ 
tion spectrum which is directly proportional to the volume 
integrated supernova rate (see YEGZ for further details). 

Accounting for the production of secondary cosmic 


© 0000 RAS, MNRAS 000, 000-000 

























































Cosmic Rays in Arp 220 3 




Energy [GeV] 

(b) Cosmic Ray Lifetimes for Arp 220 West CND 


Figure 1. The total cosmic ray lifetime (solid line) is a combination of the energy loss time (dashed line) and the energy-independent 
advective time-scale (dotted line). Energy losses for the cosmic ray protons (red line) include ionization, the Coulomb effect, and 
pion production, and energy losses for the cosmic ray electrons (blue line) include ionization, bremsstrahlung, inverse Compton, and 
synchrotron. The calorimetric nature of the nuclei means that the total cosmic ray lifetime is essentially the energy loss time as the 
advection time is negligible. 


Table 2. Best-Fitting Parameters 




East Nucleus 


West ST 

West CND 

Efficiency, rj 

5% 

10% 

20% 

5% 

10% 

5% 

10% 

Spectral Index, p 

2.1 - 2.3 

2.2 - 2.3 

2.3 

2.2 - 2.3 

2.3 

2.2 - 2.3 

2.3 

Magnetic Field Strength, B [mG] 

4.0 - 7.5 

2.0 - 2.5 

1.5 

1.5 - 2.5 

1.0 

5.25 - 8.75 

3.5 

Wind Speed, v a dv [km s —1 ] 

0-2000 

500 1400 

1600 - 2000 

0 - 2000 

0 - 300 

0 - 2000 

0 - 300 

Absorption Fraction, / a b s 

0.5 ■ 1.0 

0.6 - 1.0 

0.7 - 1.0 

0.1 - 0.7 

0.2 - 0.7 

0.1 - 0.7 

0.2 - 0.7 

Ionized Gas Density, ri[ on [cm -3 ] 

400 - 800 

400 - 600 

400 - 500 

200 

200 

2000 

2000 

Reduced Chi-Squared, % 2 

9.6 - 11.6 

10.6 - 11.5 

11.1 11.5 

9.5 - 11.5 

9.9 - 11.3 

9.5- 11.5 

9.9 - 11.3 


rays, we use our calculations of the population of ener¬ 
getic particles to predict the radio, 7 -ray, and neutrino 
spectra. For the 7 -ray spectrum, we include both leptonic 
(bremsstrahlung, inverse Compton) and hadronic (neutral 
pion decay) emission mechanisms. For the radio spectrum, 
we i ncorporate the effects of free-free emission and absorp¬ 
tion jYoast-Hull et al]|2014bl ). As in our previous models, we 
assume that the ionized gas in the nuclei acts as a foreground 
screen that some fraction (/aba) of the emitted synchrotron 
radiation passes through. When the covering fraction is low 
(/abs ~ 0 .1 — 0 . 2 ), the radio spec trum flattens at low fre¬ 
quencies (lYoast-Hull et alj|2014bl) . and when the covering 
fraction is high (/ a ba ~ 1 . 0 ), the radio spectrum turns down 
at low frequencies (YEGZ). 


As noted above, the western nucleus in Arp 220 is best 
modelled with two separate regions: an inner circumnuclear 
disc (CND) with a surrounding torus (ST). We model the 
cosmic ray populations of the two regions independently, 
treating each region as a uniform slab. However, the effects of 
absorption (free-free and 7 ^y) on the resulting radio and 7 - 
ray emission must be considered more carefully. Absorption 
occurs within each emission region, and in the case of the 
inner CND, absorption also occurs as the emitted radiation 
moves through the external, ST (see the appendix for further 
details). 


3.2 Radio Emission 


We per form y 2 tests fol l owing t he approach described in 
YEGZ iYoast-Hull et all (l2014bl ) . Comparing against ra¬ 
dio observations for each nucleus, we vary magnetic field 
strength (B), wind speed (n a dv), ionized gas density (n ion), 
and absorption fraction (/abs). While magnetic field strength 
and wind speed both directly affect the total cosmic ray life¬ 
times, the ionized gas density and the absorption fraction 
only affect the emitted radio spectrum. The free-free emis¬ 
sion and absorption coefficients are both directly propor¬ 
tional to the square of ni on , and so, the frequency at which 
the radio spectrum flattens or turns down and the amount 
of free-free emission at high frequencies both increase with 
'ft-ion* 


Observations 


Downes fe Solomon! 


iRodrfguez-Rico et al.l (120051 ) and Barcos-Munoz et al.l 
(120151) separate the integrated fluxes of the eastern and 
western nuclei from the total flux, allowing us to constrain 
parameters for each nucleus individually. As we do not 
have radio observations which are separable between the 
two regions of the western nucleus, we cannot constrain 
the magnetic field strength in each region separately. We 
therefore assume that the ratio between the magnetic 
field strength of the inner and outer regions of the western 
nucleus is equal to the square root of the ratio of t he average 
gas densities, -Bcnd/-Bst = \ J «cnd/ust ~ 3.5 (ICrutched 
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4 Yoast-Hull, Gallagher, & Zweibel 



(a) East Nucleus, 77 = 0.05 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Magnetic Field Strength ( B ) [mG] 


Magnetic Field Strength (Bond) [mG] 
2 4 6 8 10 12 14 16 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
Magnetic Field Strength ( Bst ) [mG] 

(b) West Nucleus, 77 = 0.05 


Magnetic Field Strength (Bond) [mG] 
2 4 6 8 10 12 14 16 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
Magnetic Field Strength (Bst) [mG] 


(c) East Nucleus, 77 = 0.1 (d) West Nucleus, 77 = 0.1 

Figure 2. Contours of y 2 demonstrate that changes in wind speed have little to no impact on the radio spectrum. Acceptable fits for 
no-escape (zero wind speed) cases indicate that the nuclei are proton calorimeters, although calorimetry fraction decreases with increasing 
wind speed. The best-fitting models for the radio spectrum require milligauss strength magnetic fields and are limited to a single value 
for acceleration efficiencies of 10%: 2.5 mG for the eastern nucleus and 1.0 mG for the western nucleus ST. In the eastern nucleus, for 
the lower acceleration efficiency, a larger range of magnetic field strengths (4.5 - 7.0 mG) produce models within 3cr of the best-fitting 
model and a degeneracy between magnetic field strength and wind speed is seen. 


120121) . Thus, our magnetic field strength determination for 
the innermost western nucleus is an estimate that is guided 
by Milky Way observations. 

When assuming the standard 10% cosmic ray accelera¬ 
tion efficiency, we find that agreement between the models 
and the observed radio data occurs only in a very narrow 
area of parameter space (see Fig. 2). The best-fitting mod¬ 
els for the nuclei have magnetic field strengths limited to 1.0 
mG for the western nucleus (3.5 mG in the CND, estimated 
from scaling) and 2.0 - 2.5 mG for the eastern nucleus (see 


Figs 2 and 3 and Table 2). As seen in Fig. 3, the total ra¬ 
dio emission in both nuclei flattens at low frequencies, and 
in the eastern nucleus, the radio spectrum may be turning 
over completely. This flattening of the radio spectra requires 
moderate to high absorption fractions of 50 100% in the 

eastern nucleus and low to moderate absorption fractions of 
10 - 70% in the western nucleus. 

In addition to moderate absorption fractions in each 
nucleus, we also find a high contribution from thermal emis¬ 
sion to the total radio spectrum (see Fig. 3), particularly 
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Frequency [GHz] 


(a) East Nucleus 


(b) West Nucleus 


Figure 3. Best-fitting YEGZ models for the radio spectra in the eastern and western nuclei. The total radio flux (black, solid line) in 
eastern nucleus includes non-thermal-synchrotron emission, both absorbed (blue, dotted line) and unabsorbed (red, dashed line), and 
free-free emission (green, dot—dashed line). In the western nucleus, the total radio spectru m (black, solid line) i ncorp orates emission 
from the inner CND (black, dotted line) and th e ST (black, dashed line). D ata are taken from lDownes fe Solom on (199q, black squares), 
Rodriguez-Rico et al. ( 20051. green circles), and iBarcos- Munoz et akl I201RL blue triangles). 


in the western CND where the majority of the radio emis¬ 
sion is thermal above GHz. In part, this unusually high 
fraction of thermal emission is due to the inability of the 
model to effectively fit for free-free absorption and free- 
free emission simultaneously as se en in the eastern nucleu s 
and in previous work (see YEGZ; [Yoast-Hull et al.li2014bl ). 
The ability of the models to accurately fit the fraction of 
thermal emission is further strained by the complicated na¬ 
ture of the western nucleus and the lack of separable radio 
observations. Thus, in this particular case, the fractions of 
thermal emission in the best-fitting models have limited sig- 
nificance and do not n e cessar ily contradict observations by 
lAnantharamaiah et alj (|200d ) which indicate more modest 
amounts of thermal emission (~ 10 — 20 %). 


3.3 Molecular Gas Mass 

The western nucleus Arp 220 has a very complex struc¬ 
ture which we greatly simplified. While our two-zone den¬ 
sity distribution reproduces the observed peak column den- 
sit y of ~ 10 25 cm 2 , th e mass is lower than that estimated 
bv IScoville et alj (120151 ) who derive the western nucleus gas 
mass from observations by the Atacama Large Millimeter 
Array (ALMA) of the submillimetre dust luminosity. This 
yields a total gas mass of 4.2 x 10 9 M@ for the western 
nucleus or four times our adopted value. We therefore ex¬ 
plored the effect of increased gas mass on our model by 
tripling the mass in the western torus and ran a limited 
suite models with fixed parameters. We set the spectral in¬ 
dex to p = 2.3 and ran y 2 tests over the entire range of 
magnetic field strengths and wind speeds but over a subset 
of the previously tested ionized gas densities and absorption 
fractions. We ran tests on the western nucleus for accelera¬ 
tion efficiencies of 5% and 10%. 

In comparing the results of these models with a larger 
gas mass, we find that none of the tested models are within 
3 a of the best-fitting model at the lower assumed gas mass. 
Further more, these results yield extremely short cosmic 
ray electron lifetimes such that the physical validity of the 


models are in question. I n addi tion to the higher gas mass 
estimates, IScoville et alj (120151 ) also propose a geometry 
where the molecular gas in both nuclei is confined to a thin 
(«10 pc) disc. The evolution of supernovae and cosmic ray 
interactions in this type of high molecular mass structure is 
beyond the scope of this study which is designed to estimate 
cosmic ray interaction rates in Arp 220 but will need to be 
considered when Arp 220 is detected in 7 -rays. 


3.4 Acceleration Efficiency and Proton 
Calorimetry 

In our earlier works, we demonstrated that for a given ISM, 
the YEGZ models a re highly sensitive to th e total flux of 
cosmic rays (YEGZ; lYoast-Hull et al]|2014bl ). This flux is 
primarily effected by the original energy input into cosmic 
rays and the advective time-scale, or escape fraction. The en¬ 
ergy input into cosmic rays is determined by the supernova 
rate and the assumed acceleration efficiency (rj). Within the 
uncertainty in the supernova rate, we vary acceleration effi¬ 
ciency from 5 to 20%. 

As shown above, for the standard 10% efficiency, the 
resulting best-fitting models are highly constrained in mag¬ 
netic field strength, and we find that this is also true for an 
acceleration efficiency of 20% (see Table 2). However, for a 
lower acceleration efficiency of 5%, equivalent to a lower su¬ 
pernova rate, we find a much larger range of acceptable fits 
in the eastern nucleus with magnetic field strengths rang¬ 
ing from 4 to 7.5 mG and wind speeds spanning our entire 
tested range. As such, the best-fitting models for Arp 220 
are essentially independent of wind (advection) speed (see 
Figs 1 and 2). 

In contrast, a galactic wind was a vital component in 
modelling the cosmic ray populations of the starburst galax¬ 
ies M82 and NGC 253 such that an extremely limited range 
of wind speeds resulted in fits within 3a of the best-fitting 
models. The wind speed determines the advective timescale 
for a galaxy and the fraction of cosmic rays which escape. 
Thus, wind speed is intrinsically tied to the proton calorime- 
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(a) 7 -Ray Flux, East Nucleus 


(b) 7 -Ray Flux, West Nucleus 



(c) 7 -Ray Flux, Total 



Figure 4. 7 -ray and neutrino spectra for the parameters corresponding to those for models within 3 a of the best-fitting radio models. 
The bands of spectra shown are displayed separately for each nucleus in the upper plots, where 7 = 5% is shown in red, 7 = 10% in green, 
and 7 = 20% in blue in the upper panel. In the lower plots, each summation spectra the two nuclei are displayed such that tie = 5% 
combined with r/yy = 5% is shown in red, 7 e = 5% combined with r/yy = 10% and 7 e = 10% combined with r/yy = 5% are shown in 
orange, 7 e = 10 % combined with r/yy = 10 % is shown in green, 7 = 20 % combined with rjyy = 10 % is shown in yellow, and t)e = 20 % 
combined with 7 yy = 5% is shown in purple. The black line shown on the 7 -ray plots is the differential sensitivity limit for 50 h of 
observation with the southern CTA array, taken from https://portal.cta-observatory.org/Pages/CTA-Performance.aspx. 


try fraction for a galaxy which is closely related to the total 
radio and 7 -ray emission from a galaxy. Other models for 
Arp 220 have assumed fixed advection time-scales, thus en- 
suri ng proton calorimetry with the high ga s densities in Arp 
220 (lTorresll2004l : lLacki fc Thompson! 20131 '). While our mod¬ 
els agree with others in finding that the starburst regions of 
M82 and NGC 253 are only ~40 - 60% proton calorimeters, 
we find that Arp 220’s nuclei are 65 - 100% (eastern) and 
90 - 100% (western CND) proton calorimeters (see Fig. 1). 

3.5 Future Detection of 7 -Rays and Neutrinos 

Despite the uncertainty in the calorimetry fraction and the 
total cosmic ray flux in the eastern nucleus, we can still use 
our best-fitting models to make a prediction on the emit¬ 
ted 7 -ray and neutrino fluxes from Arp 220. To calculate 
the possible 7 -ray flux, we apply the parameters of models 
within 3cr from our best-fitting radio model. Combining each 
possible set of models from the eastern and western nucleus, 
we find that the resulting 7 -ray spectra peak around ~0.3 
GeV with a maximum flux of ~ 10 10 GeV cm -2 s -1 (see 
Fig. 4). 


While this is roughly an order o f magnitude lower than 
previous upper limits for Arp 220 (ILacki et al.l 120111 ) and 
Fermi 's differential sensitivity for four years of observations, 
it is only a factor of ~2-3 times s maller than the flu x level of 
the recently detected NGC 2146 dTang et al.ll2014h . We also 
compared our 7 -ray flux with the differential sensitivity 50 
h of observations with the future southern CTA array (see 
Fig. 4) and find it to be only a factor of a few larger than 
our maximum flux. Thus, Arp 220 may still be detectable 
by Fermi within the next several years and is a good target 
for CTA, especially for energies near 1 TeV. 

In addition to making a prediction for the 7 -ray spec¬ 
trum, we can use our same results from the radio emission 
to predict the neutrino flux from Arp 220. Proton-proton 
interactions are responsible for the creation of secondary 
pions, both neutral and charged. While the neutral pions 
decay into 7 -rays, the charged pions decay into a neutrino 
and a muon which further decays into a secondary electron 
or positron and two more neutrinos. The spectrum of the 
first neutrino from the decay of the charged pion is what 
we calculate here, as the calculation of the spectra of neu¬ 
trinos produced during muon decay is more complex (see 
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Table 3. Energy Density Distribution in Galaxies 



Supernova 

Average Gas 

Cosmic Ray 

Radiation Field 

Magnetic Field 

Magnetic Field 


Power 

Density 

Energy Density 

Energy Density 

Energy Density 

Strength 


(erg yr -1 ) 

(cm -3 ) 

(eV cm -3 ) 

(eV cm -3 ) 

(eV cm -3 ) 

(**G) 

Milky Way 

2 x 10 48 

1 

1.4 

0.3 

0.9 

6 

M82 

7 x 10 48 

260 

470 

490 

2200 

300 

Arp 220 East 

7 x 10 49 

7700 

1100 

40 000 

1.0 x 10 6 

6500 

Arp 220 West CND 

1.3 x 10 5 ° 

42 000 

2500 

440 000 

1.2 x 10 6 

7000 


Notes: Values for the Milky Way are taken from Table 1.5 moraine (2011J. For M82 and Arp 220, the values for the 
cosmic ray and magnetic field energy densities taken from our best-fitting models (see Section 3 above and YEGZ). 


lYoast-Hull et al1l2014bl) . The flux of our maximum model 
is roughly 6 x 10 -12 GeV cm ~ 2 s _1 at 0.1 PeV, and at 
this energy, the range of possible models spans an order of 
magnitude in flux (see Fig. 4). Current point source sensi¬ 
tivity limits for the northern sky for IceCu be are ~ 1 0 - 9 
GeV c m ~ 2 s -1 , assuming a spectrum of E~ 2 dAartsen et al.l 
l2014bh . Thus, it seems unlikely that Arp 220 will be detected 
as a point source during a similar time frame by IceCube. 
However, extreme ULIRGs such as Arp 220 should make 
a significant contribution to a diffuse neu trino background 
dAartsen et alj|2014al : iMurase et al]|2013l 'l. 


3.6 7 —y Absorption 

In addition to accounting for 7 -ray and neutrino emission 
in Arp 220, we have also take into account the effects of 
7 —y ab sorption due to the intense radiation fields in the 
nuclei d Torres! 120041 '). At TeV energies and above, 7 -rays 
and infrared photons can interact to produce an electron / 
positron pair dDermer fe Menonll2009i ; lBottcher et al.ll2012l) . 
The resulting electrons will be of TeV energies and most of 
their energy will be lost t o emission of synchrotron X-rays 
dLacki fe Thompsonll2013h . 

Beginning at ~2 - 5 TeV, the opacity for 7 —y absorption 
in both nuclei is significantly greater than 1. This results in a 
steepening of the predicted 7 -ray spectrum at high energies 
(see Fig. 4). We find no such increase in slope in the neutrino 
flux as the steepening is an effect of interactions between the 
7 -ray and the ambient radiation field and not the cosmic 
ray proton population. Therefore, in the case of Arp 220 
and other such ULIRGs, the TeV 7 -ray flux is an unreliable 
indicator of neutrino flux. If the effects of spectral steepening 
by 7 —y absorption are accurate, then Arp 220 is unlikely 
to be detected by CTA or other ground based Cherenkov 
telescopes above ~10 TeV. 


Our maximum 7 -ray spectru m is a f actor of a 2 - 5 lower 
than previous predictions bv lTorr^ d2004 ) and less than a 
factor of 2 lower than those by Lacki fc Thompson! d2013ll . 
While the predicted 7 -ray flux will likely be detected by 
Fermi in the future, under our model assumptions Arp 220 
is unlikely to be detected as a high energy neutrino point 
source with the current IceCube observatory. Additionally, 
7 —y absorption of the TeV energy 7 -rays make the TeV 7 - 
ray flux a poor indicator of the neutrino flux in ULIRGs and 
other such systems with extremely intense infrared radiation 
fields. 

In addition, we find that milligauss strength magnetic 
fields are still necessary to reproduce the observed radio 
fluxes from the starburst nuclei, even having assumed larger 
supernova ra t es than previous models by factors of 2 - 5 
(I Torres! |2004 lLacki fe Thompsonl 120131 '). Differences in as¬ 
sumed volume and radiation field energy density across the 
various models account for the similar best-fitting magnetic 
field strengths despite the range in assumed supernova rates. 
The origins of milligauss strength magnetic fields in extreme 
starbursts and their impact on the evolution of these systems 
merit further examination. 

While the energy density in both magnetic and radia¬ 
tion fields is up from starbursts like M82 by two to three 
orders of magnitude, the change in the ratio of their energy 
densities is up by less than an order of magnitude (see Table 
3). Conversely, we see a much larger change in the ratio of 
magnetic field energy density to cosmic ray energy density. 
Because the cosmic ray energy density depends on the par¬ 
ticle energy loss rate, it does not increase at the same rate 
as the magnetic and radiation field energy densities (Yoast- 
Hull, Gallagher, Zweibel, in preparation). Thus, the mag¬ 
netic fields exceed energy equipartition with the cosmic rays 
by more than two orders of magnitude (see Table 3). 


4 DISCUSSION AND CONCLUSIONS 

In applying the YEGZ models to Arp 220, we find that the 
central starburst regions of Arp 220 are moderate to com¬ 
plete cosmic ray proton calorimeters. As such, the leptonic 
cosmic ray population is dominated by secondary electrons 
and positrons. The majori ty of these secon daries are pro¬ 
duced at low energies (e.g.. lLacki et al1l20ld ) and are likely 
a major contributor to heating of the ISM v ia ionization 
dPapadopoulosiboid : IPapadopoulos et alJhoi l). 

Based on our best-fitting models for the radio spectrum, 
we make predictions for both the 7 -ray and neutrino fluxes. 
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APPENDIX A: TWO ZONE MODELS 

Our single-zone model uses a si mple solution to the radia¬ 
tive transfer equation of (e.g., iRvbicki fe Lightmanl 1 19791 : 
iBottcher et al.ll2012l : lGhisellinill2013h 

F° bs = Fl nt x ^1, (Al) 

T 

where F^ nt is the radiative flux prior to absorption, F° bs is 
the radiative flux after absorption, and r is the optical depth 
for either free-free absorption or 7-7 absorption. This is still 
the solution for the eastern nucleus and the surrounding 
torus in the western nucleus. In the western CND, we must 
account for a standard emission and absorption region with 
an additional, external absorbing region. This observed flux 
is given by 

/ 1 _ p T inner \ 

F° bs = ( Fl nt x -— - x e ~ Touter , (A2) 

where T^ nt is still the radiative flux prior to absorption, Tinner 
is the optical depth for 7-7 or free-free absorption in the 
emission region, and r ou t er is the optical depth in the exter¬ 
nal, absorbing region. 


REFERENCES 

Aalto S., Wilner D., Spaans M., Wiedner M. C., Sakamoto 
K., Black J. H., Caldas M., 2009, A&A, 493, 481 
Aartsen M. G., et al., 2014a, Phys. Rev. Lett., 113, 101101 
Aartsen M. G., et al., 2014b, ApJ, 796, 109 
Anantharamaiah K. R., Viallefond F., Mohan N. R., Goss 
W. M., Zhao J. H., 2000, ApJ, 537, 613 
Barcos-Munoz L., et al., 2015, ApJ, 799, 10 
Bottcher M., Harris D. E., Krawczyski H., 2012, Relativis¬ 
tic Jets from Active Galactic Nuclei. Wiley, Weinheim 
Crutcher R. M., 2012, ARA&A, 50, 29 
Dermer C. D., Menon G., 2009, High Energy Radiation 
from Black Holes. Princeton Univ. Press, Princeton, NJ 
Downes D., Eckart A., 2007, A&A, 468, L57 
Downes D., Solomon P. M., 1998, ApJ, 507, 615 
Draine B. T., 2011, Physics of the Interstellar and Inter- 
galactic Medium. Princeton Univ. Press, Princeton, NJ 
Ghisellini G., 2013, Vol. 873: Radiative Processes in High 
Energy Astrophysics. Springer-Verlag, Berlin 
Lacki B. C., Thompson T. A., 2013, ApJ, 762, 29 
Lacki B. C., Thompson T. A., Quataert E., 2010b, ApJ, 
717, 1 

Lacki B. C., Thompson T. A., Quataert E., Loeb A., Wax- 
man E., 2011, ApJ, 734, 107 
Lisenfeld U., Volk H. J., Xu C., 1996, A&A, 314, 745 
Lonsdale C. J., Diamond P. J., Thrall H., Smith H. E., 
Lonsdale C. J., 2006, ApJ, 647, 185 
Mundell C. G., Ferruit P., Pedlar A., 2001, ApJ, 560, 168 
Murase K., Ahlers M., Lacki B. C., 2013, Phys. Rev. D, 88 , 
121301 

Papadopoulos P. P., 2010, ApJ, 720, 226 
Papadopoulos P. P., Isaak K., van der Werf P., 2010, ApJ, 
711, 757 


Papadopoulos P. P., Thi W.-F., Miniati F., Viti S., 2011, 
MNRAS, 414, 1705 

Parra R., Conway J. E., Diamond P. J., Thrall H., Lonsdale 
C. J., Lonsdale C. J., Smith H. E., 2007, ApJ, 659, 314 
Rodrfguez-Rico C. A., Goss W. M., Viallefond F., Zhao J.- 
H., G’omez Y., Anantharamaiah K. R., 2005, ApJ, 633, 
198 

Rovilos E., Diamond P. J., Lonsdale C. J., Smith H. E., 
Lonsdale C. J., 2005, MNRAS, 359, 827 
Rybicki G. B., Lightman A. P., 1979, Radiative Processes 
in Astrophysics. Wiley, New York 
Sakamoto K., et al., 2008, ApJ, 684, 957 
Scoville N., et al., 2015, ApJ, 800, 70 
Smith H. E., Lonsdale C. J., Lonsdale C. J., Diamond P. J., 
1998, ApJ, 493, L17 
Soifer B. T., et al., 2000, AJ, 119, 509 
Tang Q.-W., Wang X.-Y., Tam P.-H. T., 2014, ApJ, 794, 
26 

Torres D. F., 2004, ApJ, 617, 966 
Tunnard R., et al., 2015, ApJ, 800, 25 
Wilson C. D., Rangwala N., Glenn J., Maloney P. R., 
Spinoglio L., Pereira-Santaella M., 2014, ApJ, 789, L36 
Yoast-Hull T. M., Everett J. E., Gallagher J. S., Ill, 
Zweibel E. G., 2013, ApJ, 768, 53 (YEGZ) 

Yoast-Hull T. M., Gallagher J. S., Ill, Zweibel E. G., Ev¬ 
erett J. E., 2014a, ApJ, 780, 137 
Yoast-Hull T. M., Gallagher J. S., Ill, Zweibel E. G., 2014b, 
ApJ, 790, 86 


© 0000 RAS, MNRAS 000, 000-000 









